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We study the electron self-energy in a strong magnetic field when the parameter n = 
(a/2n) ln 2 (e_B/mo) ~ 1 and explore the transition between the perturbative regime n << 1 and 
the nonperturbative massless QED regime rj » 1 where chiral symmetry is broken spontaneously 
and electrons acquire the dynamically induced mass. 
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I. INTRODUCTION. 



Recently, it has been shown that, in the massless QED, the chiral symmetry of the lagrangian is broken spon- 
taneously when the system is placed in an external magnetic field ^-£§]. Spontaneous chiral symmetry breaking is 
. associated with the appearance of the gap in the electron spectrum 

m dy „(B) oc %/lBexp|-|y^| , (1) 

\o : 

the formation of the fermion condensate < ee >b~ eBm^y^B) Q, and the appearance of the massless "pion" 
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FIG. 1. 

On the other hand, in the real World with a finite physical electron mass too and when an external magnetic 
field is weak enough, we are in the perturbative regime and the shift of the pole of the electron propagator can 
be found by calculating the one-loop graph in Fig. 1 in the presence of the field. Assuming eB/m^ ^> 1, but 



(a/2n) log (eB/ml) < 1, the result is 



m(B) = toq 



eB 

47T ~~° TOn 



1 + — log 2 — + ■ 



(2) 



In this paper, we explore the region when the dimensionless parameter rj = (a/2n)\og 2 (eB/mQ) is of order 1. We 
perform an accurate summation of all leading double logarithmic contributions oc [a\og 2 (eB /m^)] and show that 
the result of the summation is a nontrivial function which has a singularity at some value of 77. This signalizes 
the breakdown of the perturbation theory and the transition to the nonperturbative regime involving spontaneous 
breaking of chiral symmetry. We solve then a self-consistent Hartree-Fock equation and find the effective mass 
dependence in the whole region of the parameter rj which interpolates smoothly between perturbative result (|^) (and 
its leading logarithmic version) at small 77 and the nonperturbative asymptotics (Q) at large rj. 



*Permanent address: ITEP, B. Cheremushkinskaya 25, Moscow 117259, Russia. 
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II. GROSS - NEVEU MODEL. 



Before going over to four dimensions, let us carry out the same program in the two-dimensional Gross-Neveu model 
The model is exactly solvable for large number of fcrmions and its dynamics is very well understood. We will 
see later that QED4 in an external magnetic field and the Gross-Neveu model have very much in common and the 
physics of the transitional region in these two theories is essentially the same. 
The lagrangian density of the model reads 



2 



4>ilpi) 2 - (2J ^il^if 



(3) 



where the fermion field carries an additional "color" index i = 1, . . .N. In the absence of bare mass mo the lagrangian 
(H) possesses U(l) chiral symmetry ipi — > e te75 ipi. In two dimensions, the coupling g is dimensionless. We will assume 
it to be small. When the number of colours N is large, the perturbative series involves only the so called daisy or, 
better to say, "cactus" graphs (see Fig. 2); all other graphs are suppressed as inverse powers of N. 



FIG. 2. A cactus graph. 

The 1-loop correction to the bare fermion mass mo is given by the graph in Fig. 3a. A trivial calculation gives 

(4) 



1 + a log — + o(a) 
mo 



where a — g 2 N/ir and A is an ultraviolet cutoff. Next, we can sum up all the leading logarithmic terms 
~ [a log(A/mo)]". They are given by the graphs in Fig. 3b. The result is 

1 — a log — 

to mo 

The leading logarithmic expression blows up at a finite cutoff which signalizes the breakdown of the perturbation 
theory even in its improved form. In this simple case, however, we can sum up not only the leading logarithmic cactus 
graphs, but just all of them if writing down the self-consistent gap equation 

, A 

m = mo + am log — , (6) 
m 

The solution of this equation is a smooth function of A. When a log(A/mo) is small, it coincides with the perturbative 
unimproved, and further, improved results (|J), (J^j. When alog(A/mo) 3> 1, we find ourselves in a nonperturbative 
regime. In this case, one can neglect the bare mass altogether, and the solution of Eq.(^) is just 
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a) b) c) 

FIG. 3. Graphs describing mass renormalization in the Gross-Neveu model in aj one loop, b) leading logarithmic approxi- 
mation, c) next-to-leading logarithmic approximation. 

to = Ae~~. (7) 

A finite gap (Q) signalizes the spontaneous breakdown of the chiral symmetry in the massless Gross-Neveu model in 
the leading order in 1 /N expansion 00. 

It is instructive to rewrite the transcendental equation (j^ in a differential form 

d y ^ (n\ 1 ra\ 

-77 = t— — > y(°) = !> ( 8 ) 

ell; I + ay 

where we introduced y = to/too and £ = logA/mo. This is nothing else as the standard renormalization group 
equation for a mass. If we now expand the right hand side of (||) keeping only the leading in a term, we come at the 
singular solution (||). 

In the next order in a we have the equation which takes care of all the leading and next-to-leading logarithmic 
terms: 

y' = ay 2 - a 2 y 3 . (9) 

In the diagram language, it corresponds to the summation of a particular infinite set of cactus graphs with one chain 
(like in Fig. 3b) or with two or more chains branched out of one and the same link where a logarithm is lost (see Fig. 
3c). The solution tends to a "fixed point" y„ — 1/a when £ — > oo. 

In the order ~ a 3 , the equation summing leading, next-to-leading, and next-to-next -to-leading logarithms reads 

y' = ay 2 - a 2 y 3 + a 3 y\ (10) 

In the graph language that corresponds to taking into account also the graphs with 2 branchings. The solution of Eq. 
( |l0| ) blows up at a finite £ like it did in the leading order. In the still next order ~ a 4 , we obtain again the equation 
with (the same !) fixed point etc. 

The solution to the exact equation (S) does not have any fixed point, neither does it blow up at a finite £. We 
plotted it (for a = .1) together with the perturbative approximations [Eq.(|^] and the solution to the equation (^)] 
and the nonperturbative asymptotics (Q) [ y &s (£,) — exp(£ — 1/gO in our variables] in Fig. 4. 



1 The Mermin-Wagner-Coleman theorem || forbids the spontaneous breakdown of continuous symmetries in two dimensions. 
The fermion dynamical mass m is an order parameter of chiral symmetry breaking only in leading order in 1/N. In the 
next-to-leading order approximation, spontaneous chiral symmetry breaking is washed out by interactions of would-be Nambu- 
Goldstone (NG) bosons. In fact, in this model we have the realization of the Berezinsky-Kosterlitz-Thouless (BKT) phase: 
while chiral symmetry is not broken in this phase, the fermions acquire dynamical mass m, and the would-be NG boson 
transforms into a BKT gapless excitation H. 




FIG. 4. Fermion mass in the Gross-Neveu model (solid line). The upper and lower dashed lines are the result of summation 
of, correspondingly, the leading and the next-to-leading perturbative logarithms and the dotted line is the nonperturbative 
asymptotics. 



III. QED IN A MAGNETIC FIELD 

To begin with, lot us spend some time to explain how the leading logarithmic 1-loop result (||) can be derived. The 
experience and insights acquired will be useful when handling higher-loop contributions. 

To calculate the graph in Fig. 1, we need to know the electron Green function in a magnetic field. The explicit 
expression has been found by Schwinger |10| . First of all, we can write 



Gb{x,v) = cxp<{ie J A^)d^ 



d^k 
W) 4 



e -W*-v)G B (k), 



(11) 



where the integral in the phase factor is done along the straight line connecting x and y. Generally, the momentum 
space fermion propagator Gsik) is given by a complicated integral over the proper time. It can be expanded in terms 
of Landau levels (the magnetic field is directed along the third axis) 



Gfl(*)=ie*p(~J)£(-l) 



n=0 



., D n (eB,k) 
fcjj — tuq — 2eBn 



(12) 



with 



D n (eB, k) = (fcii + m 



n-'- l : 2 »i„(^-)-(i + <7 1 7 a )£n-x 



where L" are the generalized Laguerre polynomials, fc|| and k± are longitudinal and transverse momentum components, 
- fc° 7 ° - fc 3 7 3 , jfc_L = /fcV + fc 2 7 2 , k\ - k 2 G - kl k\=k\ + k\. 

Fortunately, as we will see later, we do not need the full expression ( p^ ) for our purposes, but only its asymptotic 
form in the region of small momenta k <C V eB. In this region, it suffices to retain only the lowest Landau levels 
(LLL) in the spectral decomposition of the Green function which acquires then a very simple form 



G B (k) - 



k\\ + m 
fc 2 — TO 2 , 



(1 



z 7 1 7 2 ) 



(14) 



Basically, this Green's function describes free motion in the longitudinal direction of the states with a — SB/|B| = — 1 
(the factor 1 — i7 1 7 2 filters out the states with wrong sign of the spin projection a). The LLL approximation is 
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justified when both k"]_ and k 2 are much smaller than eB. Thereby, it is more consistent to set the form factor 

e ~k ± /eB (appearing due to a finite size of the Landau orbits) to 1, but we should keep in mind that all the momenta 
integrals should be cut off at kj_, k 2 < eB. 

Substituting Eqs.(|ll|) and ( |l4| ) in the Feynman integral, taking some care when handling the phase factors oc 
exp {ie J A^d^^ , and performing the Wick rotation, the 1-loop expression for the euclidean mass operator acquires 
the form PI 



MM) 



2vr 



( B 



d 2 ku 



kf, + m 2 , 



< B 



dk\ 



k\ 



( k ll 



(15) 



We assumed here that the external electron momentum p is purely longitudinal. One can show that the dependence 
of M(p 2 ,p 2 L ) on p\ is weak when p\ <C eB. 

The integral (Eq) is double logarithmic receiving a main contribution from the region m§ <C k 2 <C k\ <C eB. That 



justifies the use of the approximate LLL expression ([141). For k > VeB, the LLL approximation is not valid, but the 
corresponding contribution in the mass operator would not involve a large logarithm log(e_B/mQ). To find the mass 



shift 



|), it suffices to set p 2 — m\ in Eq.(15) and calculate the integral with double logarithmic accuracy. In the 



following, we will need, however, also the off-shell mass operator with arbitrary pjj. An important remark is that, 
if p 2 3> ?7!q, the structure of the integrand in Eq. (|l5| ) depends on the relation between fc|| and py. There are two 
distinct regions, both of them providing a double logarithmic contribution in Mi(p^). 
A) k 2 » p 2 . In that case, the integral is reduced to 



Mf(p\) 



am,Q 
~2^T 



dkj 
k 2 



arriQ 2 e B 



(16) 



B) ml < fc| < p 2 . We have 



a mo 



pi dkl 



The sum of two contributions ([!(]) and ([l7|) is 



k 2 



dk 2 
A- 



_i_ 



i:2 I 7.2 2n ™ 2 



eB 

2 



(17) 



Mi{p\) 



amo 
~2^~ 



LI - -L 2 
2 



eB 
4 



eB 

L = log— , Z = log— . 



(18) 



FIG. 5. A rainbow graph. 



2 We have chosen the Feynman gauge, but the mass denned as the position of the pole of electron propagator at zero momentum 
should not depend on the gauge and it does not jl|j2|Jl^1 . 
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Let us consider now higher loop contributions. The leading logarithmic terms come from the rainbow graphs 
depicted in Fig. 5. The n-th order contribution for the mass operator is obtained from the contribution of the 
(n-l)-th order by the virtue of the recurrent relations 



peB d 2 k f e ^ 



dk 2 



fci+(k,|-p„)2- 

This is nothing else, of course, as the n-th approximant of the integral equation 



M(ph 



to 



2tt 



•B j2 j. 



TOS 



dk\ 



^i + (k||-P||) 5 



(19) 



(20) 



[ M(p\) = mo + Mi (pjf) + •••]. 
Again, the leading logarithmic contributions in Eqs. (19), ( p0| ) come from two distinct kinematic regions: eB 3> 
^> k 2 ^> p 2 3> to,q and eB ^ k 2 ^ ^> p 2 ^>> k 2 ^> TOq. Actually, as far as the leading logarithmic terms are 
concerned, the equation (EG) can be greatly simplified: 



M(pf) = m 



a 
2^ 



log 



e_B /" dk 



- 1 /» d]^" D 



(21) 



This equation is easy to solve. One way to do it is to find the subsequent approximants 



M 



to 



M3^)=(^J mo 



— i 4 - -L 3 + -i 3 
24 6 3 



~— L 6 + — LH - -L 3 l 3 + —LI 5 
720 120 18 15 



(22) 



and to sum them over. The shift of the electron pole mass [it is defined as the solution of the dispersive equation 
m = M(m 2 ), but to the leading logarithmic order, we may just set to = M(toq) = M(L = I)] is then given by the 
series 



to = toq 



a , 2 5 
1 + — log' 



a , 2 

- l0 S ^2 



61 / a 



720 V 2tt 



log' 



eB 



(23) 



A more nice and easy way to handle Eq.(21 ) is to notice that it is equivalent to the following second order differential 
equation: 



dp 2 . 



P\\ 



dM(p 2 Y 
dp 2 . 



M(p 2 ) 



2tt pj 



(24) 



with the boundary conditions 



The solution of Eqs.(t2J),(t!5j) is 



dM(p 2 ) 



dp 2 



= 0, 



II — 



M(p\) 



pl=eB 



to 



(25) 



cos 



M(pf) =to - 



2ir to mi 



The pole mass of the electron [defined as to = M(toq)] is 



toq 



(26) 



(27) 
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The expansion of these expressions in a reproduces, of course, Eqs. 
an analogous formula obtained in [fl3f which had the form 

fa 2 eB 
m = mo ex P — log — ^ 
V 47r Too 



3) and (E3J). Note that Eq.dgTh is different from 



(28) 



The terms ~ a in Eqs.(|27|) and (|28|) coincide with each other and with the Jancovici resul t (|2| ), but the two-loop 
coefficients are already different: we have obtained 5/24 in the expansion (^) instead of 1/8 in Jl3| . Also the behaviour 
of the functions (p7|) and (^) at large values of the argument is completely different. When too tends to zero, the 
electron mass Eq.(|27|), in contrast to Eq. (p8f) , runs into singularity at some finite too signalizing the breakdown of 
perturbation theory [cf. Eq. @ !]. 

The experience with the Gross-Neveu model taught us how to handle this problem. We have to write a self- 
consistent gap equation which sums up effectively all the relevant graphs. To this end, we just have to replace the 
mass too in the integrand in Eq. ( pl| ) and, correspondingly, in the argument of the cosine function in Eq. ( p7[ ) by the 
total mass m(B). We then come to the equation^ 



to cos 




a eB 

,^-log — 
Zir ui 



TO . 



(29) 



The solution of this equation is nonsingular in the limit too — > giving a nonzero dynamical mass (|l]) of massless 
electron in a magnetic field [||. 

Like in the Gross-Neveu case, we can rewrite Eq.(p9|) in a differential renormalization-group-like form. Introducing 
the variable y = m/mo and differentiating over £ = log(e_B/TO§) we obtain 



V = 



v\Tv 



2/(0) = 1 ■ 



(30) 



We can consider the right hand side of Eq. (J30|) as an analog of /?- function for the "charge" m/m Q . The analysis is 
exactly parallel to what we have done for the Gross-Neveu model. In the lowest order in ^/a, Eq.(pOf) is reduced to 



-2/vV - 1 



(31) 



with the solution 



cos (yff log £%) ' 

which coincides with ( |27| ) and blows up at some finite too (or B). In the next order in \f(x we get 

l2a~ 



(32) 



-yVv 2 - 1 



(33) 



The right hand side of the last equation has a "fixed point" at y = yT+ n/2a to which the solution of Eq.(|33|) tends 
as \ogeB/rriQ — > oo. If keeping the term of order (^/a) 3 in the expansion of right hand side of (p0| ) the "fixed point" 
disappears while reappearing again in order {^/a) A . On the other hand, the exact equation (pOfidoes not have any 
fixed points (other than the zero field point y = 1 ) and its right hand side behaves as y/2 at large y what means 
that the total electron mass behaves as to cx V eB in a strong magnetic field. The solution to the full Hartree-Fock 
equation ([5(j]), its perturbative approximants, and the nonperturbative asymptotics (0) are plotted in Fig. 6. 



3 Note that, in contrast to Eq. (g) in the Gross-Neveu model which summed up all the perturbative graphs in the large iV 
limit and was just exact, the status of Eq. (j^j is somewhat less certain. It is written in the mean field (or Hartree-Fock) 
approximation neglecting the effects due to momentum dependence of the mass operator in the denominator of the electron 
propagator. We assume, however, that the mean field analysis is justified here as it is justified in many other physical problems 
involving gap equations (like standard superconductivity etc). Moreover, a numerical study of the nonlinear gap equation || 
confirms the formula (111) . 
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FIG. 6. Electron mass as a function of magnetic field (solid line). The dashed and dotted lines have the same meaning as in 
Fig. 4. 

The main lesson to be learned from our analysis is that, for systems involving spontaneous chiral symmetry breaking, 
the perturbative methods (even improved by summation of the leading or subleading logarithms) have an intrinsic 
barrier: one cannot go much further beyond the point alog(A/mo) ~ 1 for the Gross-Neveu model or the point 
(a /2tt) log 2 (eB/rriQ) ~ 1 for QED± in magnetic field: the perturbative approximants either blow up at some small 
enough value of mo or display an unphysical fixed point behaviour. In the region of very large A/mo for the Gross- 
Neveu model and eB/rriQ for QED, the effects due to explicit breaking of chiral symmetry in the lagrangian become 
negligible and the correct results can be obtained by solving a self-consistent mean field nonperturbative gap equation. 

In the real QED, the value of the parameter 77 ~ 1 is reached for fields of the order ~ 10 26 G. We recall that strong 
magnetic fields (B ~ 10 24 G) might have been generated during the electroweak phase transition |Q. It has been 
speculated in [|l]j^] that the character of electroweak phase transition could be affected by generation of a dynamical 
electron mass under such strong fields. On the other hand, as is seen from Fig. 6, the nonperturbative regime becomes 
prevailing over the perturbative one for values of r\ of the order 2.35 what corresponds to magnetic fields ~ 10 32 G. 
Ambj0rn and Olesen [jl5| have claimed that even larger fields, ~ 10 33 G, would be necessary at early stages of the 
Universe to explain the observed large-scale galactic magnetic fields. 

It was shown in |1| that the critical magnetic field, required for the effect of chiral symmetry breaking to be 
important, can be substantially decreased by taking into account other interactions, for example Yukawa couplings. 
It would be interesting to extend the results obtained in present paper to Standard Model. 
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